Novoa-Aponte et al. BMC Structural Biology 2012, 12:25 
http://www.biomedcentral.eom/1 472-6807/1 2/25 



^BMC 

Structural Biology 



RESEARCH ARTICLE 



Open Access 



In silico identification and characterization of the 
ion transport specificity for P-type ATPases in the 
Mycobacterium tuberculosis complex 

Lorena Novoa-Aponte 1 , Andres Leon-Torres 1 , Miyer Patiho-Ruiz 1 , Jenifer Cuesta-Bernal 1 , Luz-Mary Salazar 1 , 
David Landsman 2 , Leonardo Marino-Ramirez 2,3 " and Carlos-Yesid Soto 1 



Abstract 

Background: P-type ATPases hydrolyze ATP and release energy that is used in the transport of ions against 
electrochemical gradients across plasma membranes, making these proteins essential for cell viability. Currently, the 
distribution and function of these ion transporters in mycobacteria are poorly understood. 

Results: In this study, probabilistic profiles were constructed based on hidden Markov models to identify and 
classify P-type ATPases in the Mycobacterium tuberculosis complex (MTBC) according to the type of ion transported 
across the plasma membrane. Topology, hydrophobicity profiles and conserved motifs were analyzed to correlate 
amino acid sequences of P-type ATPases and ion transport specificity. Twelve candidate P-type ATPases annotated 
in the M. tuberculosis H37Rv proteome were identified in all members of the MTBC, and probabilistic profiles 
classified them into one of the following three groups: heavy metal cation transporters, alkaline and alkaline earth 
metal cation transporters, and the beta subunit of a prokaryotic potassium pump. Interestingly, counterparts of the 
non-catalytic beta subunits of Hydrogen/Potassium and Sodium/Potassium P-type ATPases were not found. 

Conclusions: The high content of heavy metal transporters found in the MTBC suggests that they could play an 
important role in the ability of M. tuberculosis to survive inside macrophages, where tubercle bacilli face high levels 
of toxic metals. Finally, the results obtained in this work provide a starting point for experimental studies that may 
elucidate the ion specificity of the MTBC P-type ATPases and their role in mycobacterial infections. 
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Background whereas the majority of cases (approximately 90%) 

Tuberculosis (TB) is one of the most important chal- progress to a non-infectious disease or latent TB, 

lenges in public health maintenance throughout the where mycobacteria survive in a dormant state inside 

world. According to the World Health Organization immune cells [4], Individuals with latent TB may be 

(WHO), 8.5-9.2 million new TB cases were estimated to asymptomatic during prolonged periods of time; however, 

have occurred in 2010 [1], and 1.2-1.5 million deaths TB can reactivate when the host immune response 

were caused by species of the Mycobacterium tubercu- diminishes due to malnutrition, steroid use, and HIV 

losis complex (MTBC) that includes M. tuberculosis, co-infection [5]. 

M. bovis, M. bovis BCG (vaccine strain), M. africanum, The emergence of multidrug and extensively drug- 

M. microti, M. canettii, and M. pinnipedii, which pro- resistant tuberculosis strains (MDR-TB and XDR-TB) 

duces TB in humans and some animal hosts [2,3]. Part and the lack of drugs against latent TB have become 

of the infected population will develop active TB, serious problems for TB control. Therefore, the identi- 
fication of new therapeutic targets useful in the devel- 

opment of novel drugs and vaccines against latent TB 

is essential. New anti-TB drugs, such as diarylquino- 
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essential membrane proteins that affect mycobacterial 
viability [6]. Thus, antimicrobials designed against pro- 
teins of plasma membrane are ideal because they avoid 
problems related to membrane permeability. 

Ion transport in bacteria is carried out by enzymatic 
systems that belong to either the P-type ATPase, ATP 
binding cassettes (ABC transporters) and metallic 
ion/H + -antiporter systems [7]. In general, ATPases help 
maintain the ion gradients responsible for cell volume 
control and transport of nutrients across the cell mem- 
brane [8-11]. ATPases hydrolyze ATP, releasing energy 
that is used in the transport of ions against electrochem- 
ical gradients in plasma membranes. The enzymatic 
mechanisms of P-type ATPases were initially described in 
eukaryotic cells [12-14]. These enzymes have five different 
functional and structural domains: three of these domains 
are cytoplasmic (A, actuator; N, nucleotide binding and P, 
phosphorylation), and the other two are embedded in the 
membrane (T, transport and S, class specific support 
domain) [15]. P-type ATPases have the following two 
conformational states: El, which binds ion substrates 
on one side of the cell membrane, inducing their auto- 
phosphorylation and generating a new conformational 
state; and E2, which has a lower affinity for substrates and 
therefore releases them to the other side of the cell mem- 
brane, promoting ion transport and finally recovery of the 
El state [15,16]. 

Despite P-type ATPases share the same catalytic 
mechanism based on conformational changes in their 
five structural domains, their regulation and substrate 
affinities are different [15,16]. P-type ATPases are phylo- 
genetically classified into five subfamilies (P r P v ), and 
within these subfamilies are 10 different subtypes that 
are categorized based on the transported substrate [17]. 

In this study, probabilistic profiles were constructed to 
compare and classify all P-type ATPases of the MTBC 
based on their structural features and ion transport. 
Twelve possible P-type ATPases were detected in the 
proteome of M. tuberculosis H37Rv and from other 
members of the MTBC. The high number of heavy 
metal transporters discovered in the MTBC suggests an 
important role for P-type ATPases in M. tuberculosis 
survival within macrophages. 

Methods 

Construction of hidden Markov models (HMM) 

To obtain a representative group of each phylogenetic 
subfamily, a set of 128 well-characterized P-type 
ATPases with evidence of existence at the protein level 
were retrieved from Uniprot (Swiss-Prot section) [18]. 
Each group of sequences was aligned using the Praline 
tool (http://www.ibi.vu.nl/programs/pralinewww/) with 
the BLOSUM62 matrix and Phobius transmembrane 
structure predictor, which was developed to improve the 



multiple alignments of membrane protein sequences 
[19]. The HMM package of programs [20], available in 
the Mobyle Pasteur portal (http://mobyle.pasteur.fr/cgi- 
bin/portal.py#welcome), was used to find these types of 
pumps in the MTBC proteomes. The Hmmbuild tool 
with default settings and the multiple sequence align- 
ments was used for HMM building. The default para- 
meters of the Hmmer tools use an ad hoc position-based 
sequence weighting algorithm that makes the models 
appropriate for the identification of distant members 
of the P-type ATPases family. Consensus sequences 
were generated with the Hmmemit tool. 

Search and classification of MTBC P-type ATPases 

To date, the following 10 MTBC genomes have been 
completely sequenced and assembled (NCBI): M. africanum 
GM041182 (NC 015758), M. bovis NZ2122I91 (NC 002945), 
M. bovis BCG str. Pasteur 1173P2 (NC 008769), M. bovis 
BCG str. Tokyo 172 (NC 012207), M. canettii CIPT 
140010059 (NC 015848) and five M. tuberculosis strains, 
H37Rv (NC 000962), H37Ra (NC 009525), Fll (NC 009565), 
CDC1551 (NC 002755) and KZN1435 (NC 012943); these 
proteomes were obtained from Uniprot (http://www.uniprot. 
org/). Because strains of M. microti and M. pinnipedii had 
not been sequenced, they were not included in this study. 
HMM and Hmmsearch tool were used to find P-type 
ATPases in the MTBC proteomes. 

Topology prediction 

The topology derived from predictions of transmembrane 
segments (TMS) was made with the following six pro- 
grams: TopPred (http://mobyle.pasteur.fr/cgi-bin/portal. 
py#forms::toppred), DAS (www.sbc.su.se/~miklos/DAS/), 
TMpred (www.ch.embnet.org/software/TMPRED_form. 
html), TMHMM 2.0 (www.cbs.dtu.dk/services/TMHMM), 
HMMTOP (www.enzim.hu/hmmtop/) and Phobius 
(phobius.sbc.su.se/). All programs were used with default 
settings, except TopPred, which allows the user to provide 
information about the type of organism. TMDET and 
PPM servers were used in cases where the predictions 
based on amino acid sequences did not yield reliable 
results. These tools must be fed with PDB files to identify 
TMS based on tertiary protein structure and to generate 
3D models with TMS located into hypothetical lipid 
bilayer planes. Tertiary structure models of CtpH and Ctpl 
were made with the I-TASSER tool [21] based on the 
threading strategy and were validated with the Whatif 
package of programs (http://swift.cmbi.ru.nl/servers/html/ 
index.html). 

Hydrophobicity profile construction 

The amino acid sequence of the M. tuberculosis H37Rv 
P-type ATPases and consensus sequences generated by 
Hmmemit tool were analyzed by TMHMM 2.0. 
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Hydrophobicity profiles of consensus sequences for pre- 
viously characterized P-type ATPases were used as com- 
parison patterns. 

Conserved motif analysis 

Amino acids sequences of the identified proteins were 
manually analyzed to determine the nine conserved 
motifs typical of the P-type ATPase family, as previously 
described by Thever et al [22] and others [16,23]. 

Results and discussion 

An unusually high number of cation transporter P-type 
ATPases are present in the Mycobacterium tuberculosis 
complex 

Multiple alignments of 128 reported P-type ATPase pro- 
tein sequences from a representative group of eukaryotic 
and prokaryotic cells (obtained from a curated database 
and confirmed at the protein level) allowed the identifi- 
cation of highly conserved regions within the family and 
the classification of these members according to ion 
transport. These alignments were used as the starting 
point for the construction of HMM profiles that repre- 
sented groups of ion transporters, which were also used 
to generate a consensus sequence for each group. The 
designed HMM were then used to identify P-type 
ATPases in proteomes of the different MTBC species, 
whose genome sequences have been reported. The pro- 
posed classification for the studied sequences was based 
on the obtained scores using the HMM and the 
Hmmsearch tool. 

Because the P-type ATPases are transmembrane pro- 
teins, typical alignments using BLOSUM and PAM 
matrices were not adequate. In this study, the Praline 
tool was used because it considers the differences in 
evolutionary tendencies of transmembranal and non- 
transmembranal regions. Praline applies an adequate 
matrix for each protein region based on the previous 
prediction of TMS (made with the Phobius algorithm in 
this case). As P-type ATPases contain ion-binding motifs 
within a transmembrane domain, the strategy that we 
employed produced more reliable alignments. Twelve 
hypothetical proteins with a high probability of being 
P-type ATPase transporters of metallic cations were 
identified in each of the MTBC proteomes. All proteins 
identified in the MTBC share at least 98% identity with 
their orthologs in the scanned proteomes; however, 
M. canettii and M. africanum displayed slight differ- 
ences in non-conserved regions. The pumps identified 
in the M. tuberculosis H37Rv proteome, CtpA, CtpB, 
CtpC, CtpD, CtpE, CtpF, CtpG, CtpH, Ctpl, CtpJ, CtpV 
and KdpB, were in agreement with the automatic annota- 
tion of probable cation transporter P-type ATPases in the 
M. tuberculosis H37Rv genome [24] and were taken as 
references to facilitate further analysis of results. 



A broad diversity of cations are potentially transported 
by Mycobacterium tuberculosis P-type ATPases 

Figure 1 shows the Hmmsearch scores obtained when 
the 16 HMM were used to find P-type ATPases in the 
M. tuberculosis H37Rv proteome. Those scores show the 
similarity between the identified sequences in the H37Rv 
proteome and the grouped sequences used for each 
HMM construction. The patterns thus obtained allowed 
classification of the 12 P-type ATPases into the following 
three groups: transporters of heavy metal (HM) cations, 
transporters of alkaline and alkaline earth metal (AEM) 
cations, and a group composed only of the (3 subunit of 
the prokaryotic K + transporter, KdpB. The highest simi- 
larity with the conventional P-type ATPases used in this 
study corresponded to KdpB, CtpF and CtpV, whereas 
CtpE, CtpH and Ctpl had the lowest scores. 

The HM group is composed of CtpA, CtpB, CtpC, 
CtpD, CtpG, CtpJ and CtpV, which may transport Cu 2+ , 
Cu + , Co 2+ , Ag + , Hg 2+ , Cd 2+ , Pb 2+ and Zn 2+ . Most of the 
proteins analyzed in this work correspond to the HM 
group (60%) suggesting that the active transport of heavy 
metal cations is relevant for the tubercle bacilli persist- 
ence, as it has been hypothesized for other prokaryotes 
and some unicellular eukaryotes [23]. Interestingly, evi- 
dence of toxic concentrations of intracellular heavy 
metals has been recently described in macrophages 
infected withM tuberculosis [25]. 

Alternatively, CtpE, CtpF, CtpH and Ctpl belong to 
the AEM group and may be involved in Na + , K + , Ca 2+ , 
H + and Mg 2+ transport. It is possible that some of these 
proteins correspond to Na + /K + or H + /K + ATPase 
pumps, but they only appeared to contain the a subunits 
in their tridimensional structure. Their non- catalytic- (3- 
subunit counterpart, which has been correlated with 
regulatory processes and assembly of P-type ATPases 
into the cell membrane of eukaryotic cells [10], was not 
found in any of the MTBC proteomes. 

KdpB is very different protein from the other MTBC 
ATPases. It shares 63% identity with the KdpB of E. coli 
(P03960), which corresponds to the (3 subunit of a K + 
transporting multimeric ATPase [22,26]. Despite the 
observation that KdpB has the characteristic DKTGTLT 
phosphorylation motif of this type of pump, it does not 
have a typical ion binding motif [26]. In conclusion, the 
proposed classification for P-type ATPases of the MTBC 
provides an initial approximation of their functional 
characterization. It is noteworthy that the constructed 
HMM in this work became a useful tool for the identifi- 
cation of P-type ATPases in other biological systems. 

Three different topologies can be adopted by P-type 
ATPases in Mycobacterium tuberculosis 

The six different algorithms used in the hydrophobicity 
analysis showed that all of the M. tuberculosis H37Rv 



Novoa-Aponte et al. BMC Structural Biology 2012, 12:25 
http://www.biomedcentral.eom/1 472-6807/1 2/25 



Page 4 of 12 






Novoa-Aponte et al. BMC Structural Biology 2012, 12:25 
http://www.biomedcentral.eom/1 472-6807/1 2/25 



Page 5 of 12 



(See figure on previous page.) 

Figure 1 Classification of P-type ATPases of M. tuberculosis H37Rv according to ion specificity. Amino acid sequences of M. tuberculosis 
H37Rv P-type ATPases were compared with the HMM profiles built from characterized P-type ATPases reported by Uniprot. Grouping is based on 
the scores obtained from Hmmseorch. HM group (blue bars) includes the P-type ATPase transporters of heavy metals (CtpA, CtpB, CtpC, CtpD, 
CtpG, CtpJ, CtpV). The EAM group (red bars) includes transporters of alkaline and alkaline earth metals (CtpE, CtpF, CtpH, Ctpl). The KdpB group 
(purple bars) includes the (3 subunit of K + P-type ATPase. The scores obtained from the Hmmseorch tool are shown in the bottom of each plot. 



P-type ATPases have an a-helix type TMS containing 
at least 17 amino acid residues, and the agreement of 
the topology results obtained using these prediction 
tools increases the confidence in the predicted TMS. 
Consensus transmembrane regions for each P-type 
ATPase were obtained if at least four of the algorithms 
gave similar results. Because three of the tools (TMHMM 
2.0, HMMTOP and Phobius) incorporate size and com- 
position restrictions in TMS, the results obtained using 
these algorithms are significant. 

The strategy used for topology analysis showed the 
following three different types of topology within the 
M. tuberculosis H37Rv P-type ATPases: type I, which 
corresponds to HM pumps with eight TMS (A, B and 
from 1 to 6), type II, which corresponds to AEM pumps 
with 10 TMS (from 1 to 10), and type III, which cor- 
responds to KdpB with 7 TMS (Figure 2). It was observed 
that all P-type ATPases have two cytoplasmic loops (small 
and large) that include the phosphorylation and ATP 
binding sites. The small cytoplasmic loop is located 
between TMS2 and TMS3, whereas the largest cytoplas- 
mic loop is situated between TMS4 and TMS5 (Figure 2) 
[22]. Additionally, as was expected, the N- and C-termini 
of these proteins are located in the cytosolic side, except 
in the case of KdpB, in which the N- terminus is located 
intracellularly and the C- terminus is outside the 
cytoplasm. 

Disagreement in the results was observed in CtpH and 
Ctpl analysis using the different prediction tools. Ctpl 
showed some hydrophobic amino acid sequences that 
could be considered part of TMS, but they do not fulfill all 
of the characteristics associated with TMS; meanwhile, 
CtpH did not show the expected hydrophobicity pattern 
for P-type ATPases (Figure 3). Therefore, at first it was diffi- 
cult to determine whether CtpH and Ctpl were in fact 
typical of P-type ATPases. Recent studies have classified 
CtpH and Ctpl proteins of M. bovis as FUPA 24 (TC No. 
3.A.3.24), i.e., "Functionally uncharacterized P-type ATPase 
family 24", and classify them as transporter proteins in the 
TCDB-Transporter Classification Database (http://www. 
tcdb.org). FUPA 24 proteins are homologous to P-type 
ATPases but have an unusually large N-terminal segment 
that makes them twofold longer than typical P-type 
ATPases [23] with two TMS. In addition, functional motifs 
of P-type ATPases are located within the C-terminal region 
of FUPA 24 [23]. 



By contrast, topology analysis in this study showed 
that M. tuberculosis H37Rv CtpH and Ctpl contain more 
than the expected TMS (from 3 to 12) for FUPA 24 pro- 
teins. To overcome the discrepancy between TCDB and 
the six topology tools used in this work, the TMDET and 
PPM servers were used. These tools allow for the identi- 
fication of TMS based on the tertiary structure of the 
protein; thus, not only can the correct size and compos- 
ition of a-helixes be guaranteed, but also their adequate 
disposition and organization in the plasma membrane 
can be determined. Ten TMS were found in CtpH and 
Ctpl using both the TMDET and PPM servers. These 
segments are located with high probability in a hypo- 
thetical lipid bilayer as shown in Figure 4. This result 
strongly suggests that CtpH and Ctpl have a type II top- 
ology, as has been determined for other members of the 
AEM group. 

Hydrophobicity profiles suggest insights into ion 
transport specificity 

The "topological clustering" based on the hydrophobicity 
profiles obtained with the TMHMM 2.0 algorithm can 
be used as a predictive tool to determine substrate speci- 
ficity for HM pumps [27]. Here, this strategy was applied 
to analyze HM, AEM and KdpB groups of M. tubercu- 
losis H37Rv P-type ATPases. The hydrophobicity pro- 
files from each M. tuberculosis P-type ATPase was 
compared with the obtained profiles from previously 
characterized P-type ATPases. As shown in Figure 3, the 
hydrophobicity profiles of CtpA, CtpB and CtpV were 
similar to the consensus Cu + P-type ATPases; the 
CtpC profile was similar to the different Zn 2+ trans- 
porters Zn 2+ , Cd 2+ /Zn 2+ and Pb 2+ /Cd 2+ /Zn 2+ /Hg 2+ . In 
addition, CtpD, CtpG and CtpJ had regions similar to 
the Cd 2+ /Zn 2+ /Co 2+ transporter. In the case of the 
AEM group, CtpF possesses TMS that were similar to 
the hydrophobic profiles obtained for consensus of Na + , 
Mg 2+ and Na + /K + ATPases. Additionally, CtpE was simi- 
lar to Na + or H + transporters; however, its hydrophobi- 
city profile is closer to that of Na + P-type ATPase. 
Considering the C-terminus region of CtpH (819 amino 
acid residues), this pump had a hydrophobicity profile 
very similar to the Ca 2+ PMCA consensus. Finally, KdpB 
had the exact same hydrophobic profile as that which 
was determined for E. coli KdpB. 
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Figure 2 Predicted topology for P-type ATPases of MTBC. Type I topology (HM group ATPases) (a), type II topology (EAM group ATPases) 
(b), type III topology (KdpB group ATPases) (c), Ion binding motifs (number 4) are indicated in the corresponding transmembrane domain. 



In general, there are the following two types of hydropho- 
bicity profiles: the first comprises the AEM and the second 
that corresponds to HM. Interestingly, the relative position 
of TMS along amino acid sequences is highly conserved be- 
tween both types of hydrophobicity profiles. The profile for 
AEM is characterized by a single TMS in the C-terminal 
region and two in the N-terminal region, whereas profiles 
for HM P-type ATPases have three highly hydrophobic 
regions, one of them in the N-terminus, one in the middle 
of the amino acid sequence, and the third in the C-terminal 



end. Although the hydrophobicity profile of consensus 
sequences is similar among enzymes in each group, there 
are significant differences between profiles that allow their 
differentiation according to the transported ion. 

P-type ATPases of the Mycobacterium tuberculosis 
complex possess the characteristic motifs involved in 
cation transport 

The common catalytic mechanism of P-type ATPases is 
partially supported by conserved core sequences that are 
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Figure 3 Hydrophobicity profiles of M. tuberculosis H37Rv 
P-type ATPases. Hydrophobicity profiles allow grouping of P-type 
ATPases with seven different ion specificities. 



useful in the recognition of these pumps in proteomes 
[15]. To find the conserved motifs for P-type ATPases 
along the protein sequences, manual analysis was carried 
out because the servers commonly used for analysis 
of conserved motifs do not identify subtle variations 
within particular motifs of P-type ATPases. Nine 
conserved motifs typical for P-type ATPases were 
found in the MTBC pumps (Figure 5). We observed 
that each motif sequence was identical within ortho- 
logs of the MTBC. From N-terminus to the C- 
terminus, the motifs were as follows: motifs 1 
[(PVA)G(DE)] and 2 [P(AS)D], related to conform- 
ational changes, together with motif 3 [TGE(SA)], 
associated with phosphatase activity, were located 
between TMS 2 and 3 in the actuator domain; the 
sequence of motif 4, which determines ion specifi- 
city, can be [PEG(LM)] or [(CSA)PCA(LV)], which 
was located at the end of the fourth TMS; the phos- 
phorylation site [DKTGTLT] was found in the P do- 
main; the remaining motifs, motif 6 [(KI)GA(PVA) 
(EDA)], an ATP binding facilitator, motifs 7 [(DV) 
(ASPI)(VP)(KAR)] and 8 [(MLV)I(TS)GD], involved 
in phosphorylation catalysts, and motif 9 [(VTC)AM 
(TV)GDG(VSAT)ND(AV)(PAL)A(LI)(RKA)(QMAD)A 
(DNT)(VI)G(VI)(AG)(MV)], the hinge motif, which 
provides the flexibility necessary to achieve conform- 
ational changes during the pumping process, could 
be found between the fourth and fifth TMS [22]. 
Motifs 8 and 9 contain amino acid residues [TGDN 
and GDGXND] responsible for the coordination of 
Mg 2+ as a cofactor of the enzymes [15,16,22,23]. 

It was observed that characteristic motifs of P-type 
ATPases exhibit slight variations compared with previ- 
ously reported sequences for eukaryotic organisms 
[22], The most conserved motif (motif 5) was almost 
the same within the 12 P-type ATPases, except CtpH 
and KdpB, which contained the conservative substitu- 
tions [DKTGTL(TS)] and [DKTGT(LI)T], respectively. 
The sequence of motif 4 [PEGL] has been previously 
associated with binding of AEM cations [22], and it 
has been found in CtpE, CtpF, and Ctpl. Variation in 
the fourth position of this motif [PEGM] in CtpH has 
not been reported previously. CtpA, CtpB, CtpC, 
CtpD, CtpJ, CtpG and CtpV have the motif 4 sequence 
[(CSA)PCA(LV)] characteristic of HM transporters. Al- 
ternatively, KdpB, which corresponds to a (3 subunit of 
a multimeric P-type ATPase, did not have motif 4. 
This finding can be explained by the observation that this 
subunit mediates phosphorylation/dephosphorylation and 
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Figure 4 The CtpH and Ctpl mycobacterial P-type ATPases (type II topology) exhibit ten TMS in their C-terminal half (different color 
helixes). These tertiary structure models were generated with the l-TASSER tool and were modified with the PPM server [89] to include the 
possible location of the lipid bilayer. Dummies (green color) correspond to the location of the carbonyl groups in the bilayer. (a) CtpH model 
generated with the PPM server, (b) Ctpl model built based on the TMDET results. 

V J 



energy transduction during K + transport. In E. coli, the co- 
ordination function lies mainly in the KdpA sub unit [26]. 

Cu + , Zn 2+ and Co 2+ pumping is mediated by P-type 
ATPases in the M. tuberculosis complex 

The P IB phylogenetic group of P-type ATPases is com- 
posed of HM pumps, and according to ion specificity, 
the P IB group is subdivided into five subgroups (1 to 5). 
Alignments of consensus sequences from characterized 
P ffi ATPases with the HM P-type ATPases of MTBC 
were used for locating additional transmembrane motifs 
of this type of transporter [28-30], as shown in Figure 6. 
CtpA, CtpB and CtpV have the characteristic motifs of 
Pib-i group members (Table 1), corroborating their ion 
transport specificity as predicted by hydrophobicity pro- 
files, and suggesting that Cu + is transported preferen- 
tially to Cu 2+ by these pumps. 

This result is expected, because the special reducing 
conditions inside mycobacteria [27] are similar to the 
intra-phagosomal environment in macrophages where 
the bacilli reside. Recently, it was reported that the 
M. tuberculosis H37Rv CtpV is a Cu + exporter P-type 
ATPase, reinforcing our predictions [31]. CtpD and CtpJ 
exhibit the [SCP] and [HEGT] motifs of the P IB _ 4 group 
in the TMS6 and TMS8, respectively. Some reports have 
indicated that amino acids involved in Co 2+ transport 
are still not well understood; for this reason, the absence 
of residue N in TMS7 cannot rule out the possibility 
that CtpD and CtpJ are Co 2+ transporters. CtpC only 
exhibits the [CPC(X) 4 S] motif in TMS6, which allows its 
classification as Pi B _ 2 . in fact, it was recently reported as a 
Zn 2+ P-type ATPase in M tuberculosis H37Rv [25]; defi- 
ciency in this pump produced zinc accumulation within the 
mycobacterial cytoplasm, resulting in impaired intracellular 



growth of tubercle bacilli [25]. Finally, CtpG does not have 
any of the additional motifs present in TMS6, TMS7 and 
TMS8; however, it possesses the [WI(YE)(RG)] sequence 
just before TMS6, between positions 406 and 409, and the 
[LS] motif located on TMS7, which is associated with Zn 2+ 
P-type ATPases [32,33]. 

All of the results from this work provide predictive 
evidence for experimental studies to establish the ion 
specificity of MTBC P-type ATPases and their role in 
mycobacterial infection. It has been observed that some 
P-type ATPases of tubercle bacilli are over-expressed 
under conditions that mycobacteria face during infection 
[25,34,35]. For example, the expression level of CtpF and 
CtpC increase when M. tuberculosis is exposed to isoxyl, 
tetrahydrolipstatine and SRI#9190 antimicrobial com- 
pounds, indicating that these pumps might contribute to 
intrinsic resistance of mycobacteria to antimicrobial 
drugs [36]. Two additional studies indicate that hypoxia 
induces upregulation of the AEM transporter CtpF 
[37,38]; this observation is in agreement with those of 
another study that reported that CtpF is part of the reg- 
ulon of the DosRS system, a relevant regulator of latency 
in M. tuberculosis [39] . In addition, upregulation of CtpF 
is observed in vitro when M. tuberculosis is incubated in 
the presence of S-nitrogluthatione GSNO, ethanol, H 2 0 2 
and nitric oxide [40]. Moreover, CtpA, CtpC, CtpG, 
CtpV and CtpF are also induced when M. tuberculosis is 
phagocytized by macrophages [25,34,35]. 

Conclusion 

Mycobacteria are unicellular organisms that respond to 
environmental stimuli, and the transport of substances 
across the plasma membrane could play a fundamental 
role in their adaptability. Computational analysis shows 
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Figure 5 Conserved motifs for the P-type ATPases identified in M. tuberculosis H37Rv. The conserved motifs were compared with the 
motifs characteristic of the P-type ATPase superfamily. Identical residues, conserved substitutions, semi-conserved substitutions and unrelated 
residues are indicated in red, blue, green and black, respectively. 
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Figure 6 Alignments of consensus sequences from characterized P| B ATPases, with the HM P-type ATPases from M. tuberculosis H37Rv. 

(a) P| B -i motifs, (b) P| B _ 2 motifs and (c) P| B - 4 motifs. 
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Table 1 Additional conserved motifs observed in heavy metal transporter P-type ATPases 
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The additional motifs are located on TMS6, TMS7 and TMS8. *, pumps without a complete set of motifs (see discussion). 



that each MTBC species has a consistent aggrupation of 
the 12 P-type ATPases involved in ion transport. In this 
context, M. tuberculosis strains H37Ra and H37Rv share 
identical sequences for P-type ATPases, facilitating sub- 
sequent genetic studies using the attenuated strain 
H37Ra. The large number of HM P-type ATPases 
expressed by the MTBC strongly suggests that they 
could be essential for the bacteria to counteract the 
increased level of HM accumulated by macrophages 
after infection with tubercle bacilli. Thus, compensatory 
ion transport strategies could be used by mycobacteria 
to survive in host cells. 

The different bioinformatics approaches used in this 
work to analyze the P-type ATPases identified in the 
MTBC are in agreement with the initial classification 
from the HMM search. The results obtained show that 
M. tuberculosis has the following three groups of P-type 
ATPases: HM transporters (CtpA, CtpB, CtpC, CtpD, 
CtpG, CtpJ and CtpV), AEM transporters (CtpE, CtpF, 
CtpH, and Ctpl) and the KdpB protein, which corre- 
sponds to the |3 subunit of a multimeric K + ATPase 
transporter exclusive to prokaryotes. Hydrophobicity 
analysis identified a-helix type TMS grouped into the 
following three topological types: type I (HM group), 
type II (AEM group) and type III (KdpB group). Interest- 
ingly, we report a possible mis-annotation for CtpH and 
Ctpl in the TCDB Database, where they are classified as 
FUPA 24 type with two TMS, unlike the ten TMS iden- 
tified for these unusually large transporters in this work. 
Finally, a counterpart of non-catalytic (3 subunits of 
Na + /K + or H + /K + ATPases does not exist within the 
MTBC proteomes. 
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